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The effects of fluid inertia and turbulence on the force coefficients of 
squeeze film dampers are investigated analytically. Both the convective and the 
temporal terms are included in the analysis of inertia effects. The analysis of 
turbulence is based on friction coefficients currently found in the literature for 
Poiseuille flow. 

The effect of fluid inertia on the magnitude of the radial direct inertia 
coefficient (i.e. to produce an apparent "added mass" at small eccentricity ratios, 
due to the temporal terms) is found to be completely reversed at large eccentricity 
ratios. The reversal is due entirely to the inclusion of the convective inertia 
terms in the analysis. 

Turbulence is found to produce a large effect on the direct damping coefficient 
at high eccentricity ratios. For the long or sealed squeeze film damper at high 
eccentricity ratios, the damping prediction with turbulence included is an order of 
magnitude higher than the laminar solution. 


NOMENCLATURE 
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C , C 
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C , C 
rt tt 


D , D 
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inner cylinder radius 

dimensionless journal center radial acceleration 
outer cylinder radius 

dimensionless damping coefficients in (r,t) directions due to tangential 
velocity 

damping coefficients = (C , C ) * pkL/6^ 

dimensionless_inertia coefficients in (r,t) direction due to normal 
acceleration a 


r 

D rr > D tr = inertia coefficients (D rr > .D ) * pkL/<$ 3 u) 
e - circular centered orbit radius 

f r » f^ = dimensionless fluid film force in (r,t) direction 
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V.*s 


= inertial wall shear stress functions 
H = 1 + ec0: dimensionless film thickness 

h = b6H = film thickness 

I 00 , l 0 ^j = momentum integrals over the film thickness 


k 


fc 0 J \ 


L 

P 


= geometry parameter = (L/b) for short SFD assumption, 1 for others 
= parameters depending on the nature of the flow 

= squeeze film damper length 
= pressure 


P = p6 /(km y) = dimensionless pressure 

q q ^ = dimensionless local flow rates in (0,g) direction 

2 2 

Re = w 6 b /v ~ squeeze Reynolds number 

h 


Rep=|^ [(q 0 +H) 2 + (L/b q ? ) 2 ] 


u 


u* 


u 


m 


-/n 
J 0 


w 


w 


w 


“ Poiseuille flow Reynolds number 
55 fluid relative velocity along lubricant film 
= absolute fluid velocity along lubricant film = u+bo) 

= u/brn: dimensionless fluid velocity 

udn= mean fluid velocity along lubricant film 
= fluid velocity across lubricant film 
= v/bdco = dimensionless fluid velocity 
= e = dimensionless journal center tangential velocity 
= fluid velocity in the axial direction 
= w/Lo) = dimensionless axial velocity 

wdn = mean fluid velocity in axial direction 
w = frequency of damper motion 

t “ time 

(x,y,z) = moving coordinate system 

(x,y,z)* = fixed coordinate system 


m 


Jo 
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(x^ 9 a^ 9 a.^ = coefficients eventually depending on Re for turbulent motion 

y = 3H/90 = film thickness gradient along circumferential direction 

— 2 

T = r (urn ) - dimensionless inertia function for long SFD assumption 
6 = (b-a) /b = clearance ratio 

e - e/bd = dimensionless circular orbit radius, eccentricity radius 
(0,tv,£) = dimensionless coordinates = (x/b, y/h, z/L) 

At ^ , At ^ = wall shear stress difference in (9,5) direction 

P = fluid density 

]i - fluid viscosity 

v - u/P - kinematic viscosity 

x = t to = dimensionless time 

Subscripts: 

0 = inertialess or purely viscous 

1 = inertial 


INTRODUCTION 


Squeeze film dampers (SFD) are designed to have a stabilizing effect on the 
rotordynamics of turbomachinery. This has generally been accomplished by using the 
Reynolds effect in a thin oil film around a bearing to produce a predictable 
damping coefficient. The increase in size and speed of modern turbomachinery using 
light viscosity oils has brought the need to include fluid inertia effects in the 
design analysis. Sparked by the recent pioneering work of Tichy [1-4], researchers 
are now extending the lubrication theory into the range where the Reynolds (slow 
flow) assumption is no longer applicable. 

At least for some simple geometries and motions, the fluid inertia effects 
have been shown to be quite significant. 

To the rotordynamicist fluid film forces and dynamic coefficients are more 
important than velocity or pressure fields. Analytical [2,7] and numerical [5,6] 
approaches have been developed for calculating the damping and inertia coefficients, 
assuming motions of small amplitude a about an equilibrium point. In this case, 
it can be shown that the convective jnertial terms may be neglected in the equations 
of motion, since they are of order a while the temporal terms are of order a. In 
all these analyses, the trend of the damping and inertia coefficients is to increase 
as the static eccentricity ratio increases, a fact that has been shown to be true 
in practice. 
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However, for large excursions of the journal center about its centered position 
the full inertial term should be retained in the momentum equations. The temptation 
to neglect convective terms in order to simplify and linearize the problem is no 
longer justifiable even for very simple cases such as the long or short bearing 
solutions. For example, reference [9] recently presented numerical calculations for 
the dynamic coefficients of an SFD performing circular orbits (CCO) about the center 
of the bearing housing. Using the same approach as in [5], the convective inertia 
effects were neglected, so the coefficients have the characteristic form described 
above. This behavior of the fluid film forces will be shown to be in error even for 
moderate eccentricities and totally incorrect at large orbit amplitudes. Furthermore, 
in reference [9] the direct damping and inertia coefficients for the cavitated SFD 
were found to be one half the value of the full film case and independent of the 
inertia parameters of the fluid. This appears unreasonable since if cavitation is 
allowed, the extent of the region where the film is broken will be influenced by the 
magnitude of the inertial forces. Our purpose in the present analysis will be to 
determine the dynamic coefficients taking into account the full inertial terms for 
simple geometries in order to understand better the action of viscous and inertial 
forces in an SFD. 

The inclusion of inertia complicates the problem in a SFD, and turbulence 
effects make the problem even more involved. Unlike the journal bearing case where 
a considerable amount of analytical and experimental work has been done, turbulence 
in squeeze film dampers remains rather obscure due to the lack of experimental data 
or a good understanding of the mechanics of squeezing flows. Nelson [11] used the 
empirical friction coefficient for pure Poiseuille flows in an attempt to include 
turbulent effects for the long SFD case. No satisfactory results were obtained 
since the fluid apparent viscocity was used to calculate the empirical friction 
factor for the flow. 

Tichy [4] suggests that turbulent flow in SFD’s will occur at higher Reynolds 
numbers than for Poiseuille flows, i.e. Rep > 2000. This assertion seems reasonable 
since the velocity field in a SFD is constantly changing and adjusts itself to the 
normal motion of the boundary, thus making the flow more stable. It also seems 
reasonable that transition from laminar to turbulent regions should be smooth in 
order to satisfy continuity of the flow. All these considerations make the problem 
more untract able and point out the urgent need of experimental data. In the mean- 
time, it will prove helpful to use the empirical correlations currently found in 
the literature and thus obtain upper bounds for the forces and dynamic coefficients 
when turbulence is present in the flow. 


STATEMENT OF THE PROBLEM 


Figure 1 shows the geometry of the SFD system. The equations of motion for 
the flow in the annular region between a whirling nonrotating inner cylinder and 
its bearing housing are stated in a moving coordinate frame. Appendix A contains 
the details of the integration of the motion equations across the lubricant film, 
to finally obtain in dimensionless form: 


Re 


— I + — i 

80 90 8? 0£ 


a -Hk -J- a t 

80 ' 0Tj 


( 1 ) 
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| 9 9 

Re le I e5 + 1? 1 kk 


- -» «>/ l > 2 H + ^ 


( 2 ) 


80 q 0 


+ 3C q C 


= 0 


( 3 ) 


This system of equations must be solved with appropriate boundary conditions 
for the flow rates (q , q ) and the pressure p. Analytical solutions to the problem 
are extremely difficult since the exact form of the wall shear stress difference 
(Atq^, At^) is unknown and some assumption regarding their functional form becomes 
necessary. Presumably, the problem may be solved numerically on a computer with its 
full tridimensional complexity, but such an effort may prove to be unnecessarily 
costly or even impractical. 


As a first approximation to a practical solution of the problem, we assume that 
for the laminar region, and even in the presence of turbulent effects, the shear 
stresses at the walls may be written as: 


Nn “ -k e f e 

H 


(4a) 


At 


?n 




Re 


(4b) 


The approximate form of the functions k Q , k r , f A and will be discussed later 

U c u c, 

in the analysis. Note also that in equations (4J we have included an explicit 
contribution of inertia to the wall shear stress difference. 


Once a solution to the system of equations (!) to (3) has been obtained the 
fluid film forces acting on the inner cylinder are calculated by integration of 
the pressure distribution over the flow region. For rotordynamics applications, 
the forces are expressed in terms of damping and inertia coefficients. Let (f^,f^) 
be the radial (along the centerline of both cylinders) and tangential dimensionless 
fluid film forces, and given by: 


f 

r 



cos 0 dR = - C! V 
rt t 

sin 0 dR = - C V t 



- D 


tr 


R = { (K 0 < 27 F , 0 <£< 1 } 


(5a) 

(5b) 


In equations (5), V t , a r are the dimensionless journal, center, tangential 
velocity and radial acceleration, respectively ; and (C^ t , C^. t ) , (D rr , D tr ) are t ^ ie 
dimensionless damping and inertia coefficients. The dimensional counterparts of 
these coefficients are given by the relations: 


C. . 
13 
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s 3 13 


D. . 
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Our interest is to obtain approximate solutions to the uncavitated case, and 
present the dynamic coefficients with inertia and turbulence effects accounted for 
in the flow. This is a necessary step preliminary to any more refined analysis, 
since it will contribute to a better understanding of the problem. We will treat 
the laminar and turbulent solutions for the long and short SFD’s separately. No 
cavitation is considered in the flow region. This last assumption is unrealistic 
for some cases, but permits us to do a first treatment of the fluid film forces 
with inertia and turbulence included, and will be accurate for high supply pressures. 
Furthermore, the direct effects of fluid inertia and turbulence will be isolated 
from the indirect effect caused by changes in the region of cavitation. 


LAMINAR FLOW SOLUTIONS 


Long Bearing Assumption 

In this section we assume that the cylinders are infinite in extent, or that 
very tight end seals are placed at the ends of the SFD, or that the axial flow q 
is negligible. We are left with the equation: ^ 


q 


0 


u 

m 




(q 0 +H) 


3T = ik [ ~ 7 + 35a i 




(7) 

( 8 ) 

(9a) 


P 


+ Re p . 

l 


Here we have divided the pressure into two parts so that p. contains explicitly the 
influence of inertia. 


For small Reynolds numbers, Re<l, Brindley [7] and the first author [10] have 
found that k = 12 as in the usual lubrication approximation, and that the inertial 
pressure gradient is given by: 




[-6 + 6 u 
m 


+ 54 u 2 ] 
m 


(9b) 


For large Reynolds numbers, assuming that the flow remains stable and laminar, 
the inviscid pressure gradient is given as: 



90 H 



(9c) 


Equation (9c) is different from the result presented by Tichy in [4], 
apparently due to an error in the boundary conditions used for the inviscid flow 
region. 
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For both extreme values of the inertial parameter Re, the dimensionless flow 


rate q is unaffected by inertia and is equal to: 
o 


l e 


= 2(e -1) 
(2+e 2 ) 


( 10 ) 


In order to make a quantitative comparison of equations (9a-9c) we let: 


H 9p i 

r =w = r V> 


(ID 


be a function that indicates the magnitude of the inertia term in the pressure 
equation. 

Figure 2 shows equations (9a-9c); curve A represents (9a) with f^ = 0, a^-1.2, 
curves B & C, equation (9b-c) , respectively. 

A brief look at figure 2 shows the surprisingly similar behavior of the flow 
for the large range of Re considered. The actual value of T for moderate Reynolds 
numbers will lie between curves A and B; curve F shows the best fitting line 
between the A and B curves from which we select the inertial contribution to the 
wall shear stress difference as: 

2 


f e = 


35 


T 


H' 


( 12 ) 


Thus, we assume that for moderate Reynolds numbers, 
difference is approximately given by: 


At 


en 


-12 


(q 0 +H) 


+ Re 


4 

35 



the wall shear stress 


(13) 


With these considerations, the pressure field for the flow can be determined, 
and from this, the dynamic coefficients. As previously stated, only the uncavitated 
SFD is treated here so that the effect of inertia on the fluid film forces can be 
clearly isolated. Otherwise, the extent of the cavitated region is dependent on 
the Reynolds number and the dynamic coefficients must be determined numerically for 
each change in the inertia parameter. 

Integration of equations (8) and (9), subject to the continuity condition for 
the pressure field, is relatively easy and is given in [8]. The dynamic coefficient 
for the long SFD with laminar flow come to be: 


C = D =0 
rt tr 


24ir 


fct (2+e 2 )(l-e 2 ) 1 ' 2 


(14) 


Note that the direct damping coefficient C is the same as in the inertialess 

solution. Analytical expressions for the direct inertia coefficient D are given 

in Table I. Figure 3 shows a comparison of this coefficient for the Hxfferent 

cases considered. The behavior of the inertia coefficient D divided by Re is 

rr 

surprisingly similar for the large range of Reynolds number considered. The largest 
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difference occurs at small and large eccentricities but is never more than 20%, 
(between Re=0 and Re=°° ). The rapid decrease of the inertia coefficient as the 
eccentricity ratio grows larger appears to contradict the recent results presented 
in [9] . The reason for the discrepancy is that in the MTI study only temporal 
acceleration effects were accounted for in the equation of motion while here both 
the convective and temporal inertia terms are retained. The upper dashed curve in 
figure 3 shows the inertia coefficient when only the temporal effects are 
included. The coefficient grows rapidly and even goes to infinity as the orbit 
size approaches the radial clearance. 


Short Bearing Assumption 

In this section we assume that the SFD has small L/D ratios, the ends are 
open to the atmosphere, and for simplicity we also assume that no high externally 
induced axial flow is present in the damper. The reason for the latter assumption 
is to avoid pressure boundary conditions which would require the explicit presence of 
the inertial parameter Re. As is current practice for the short journal bearing 
analysis, the circumferential flow is assumed to be negligible. We set k=(L/b) 
in eqs. (1) to (3) and get the following set of equations: 


q e = - H 

(15) 

35 = Y 

(16) 

H !Z° = _ 12 1l 

35 h 2 

(17) 

H _ f + !!i.!3L (q 2) 

35 5 30 H 35 w 5 ' 

(18a) 


The axial inertial pressure gradient obtained for small Re using a regular 
perturbation solution in Re is given by Tichy (3) as: 


H BP. 

w 


= 1.2 


Bq 


30 


51 B 2 
35H 3-5 vq £ J 


(18b) 


Assuming the flow remains stable and laminar, for large Re the inviscid 
pressure gradient is given by: 


H 3p. 

bT 


Bq 

*30 


1 , n 2 


(18c) 


From a quantitative comparison of equation (1.8a-c) , for moderate Reynolds 
numbers we select the inertial contribution to the wall shear stress difference 
to be: 



2 3 __ 

35H 3g 


2 1 3q 
(q 5 } + 10 36 


(19) 
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Thus, we assume that the wall shear stress difference for 
approximation is given by: 




2 3 /• 2 * 1 8q 5 

35H h (q 5 } 10 ¥T 


the short SFD 


(20) 


With these considerations, equations (17) and (18a-c) are integrated to obtain 
the pressure field. For the full film assumption, the dynamic coefficients come to 
be: 


C - D = 0 
rt tr 


tt 


(1-e 2 ) 


3/2 


( 21 ) 


Note that the direct damping coefficient C is the same as in the inertialess 
solution. Analytical expressions for the direct 1 ' inertia coefficient D are 
given in Table II. For the different cases considered, figure 4 shows r Ehe inertia 
coefficient divided by Re as a function of the eccentricity. As in the long 
SFD case, the form of the inertia coefficient is surprisingly similar for the large 
range of squeeze Reynolds numbers considered. Note the tremendous influence that 
the convective inertial terms have on the coefficient when compared to the dashed 
curve which is based only in the inclusion of temporal effects on the equation of 
motion. Thus , analyses based on small perturbation about an equilibrium point are 
in large error compared to the exact solution, if the orbit is large. 


TURBULENT FLOW SOLUTIONS 


The inclusion of turbulence effects into the flow complicates the problem 
enormously. Although the mechanism of turbulence for fully developed Couette and 
Poiseuille flows has been studied extensively, both analytically and experimentally, 
and many contributions to the analysis of flow in narrow channels have been given in 
the past years; the mechanics of squeezing flows are far more complicated. The 
subject still remains obscure due to the complete absence of theoretical-empirical 
formulation and the lack of experimental results. 

Undaunted, we assume that the coefficients & k to be used in turbulent flow 
in a SFD are given by: ** 

k r = k Q = 12 +0.005 Rep (22) 

where 

% 

Rep -|S. [(q e+H ) 2 + (L/b) 2 q ? 2 ] (23) 

is the Poiseuille Reynolds number currently found in the literature. 

Relation (22) was obtained as the best fitting curve between the experimental 
correlation given by Hirs (12) and the analytical results based on the mixing length 
theory given by Elrod and Ng (13). Here we have assumed that the transition from 
the laminar to turbulent regions in a SFD must be smooth in order to insure 
continuity of the flow. 
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The assumed expressions may be far away from the actual expressions which 
should be obtained from experimental results. However, we have chosen them in the 
absence of better empirical formulations, and the results obtained will prove to be 
upper bounds of the actual forces and dynamic coefficients. 


Long Journal-Bearing Assumption 

As in section 3. a, the axial flow is neglected and the pressure gradient field 
is given by: 

, 2 

rln Tie* (H+q~) _ q. 

H - -(12 + 0.005 f - 1H + q 0 l) j 2 - + |f * [-7+46 -f-] (24) 

* H H 

q 0 = q 0 (e. Re, 6) (25) 

Note that we have assumed that f 0 given in (12) prevails even in the turbulent 
regime. A very simple computer code was written to obtain the pressure field. Using 
numerical integration, the dynamic coefficients were calculated for an uncavitated 
SFD with a clearance ratio 6 = 0.001. 

Figure 5 shows the direct damping coefficient C as a function of the 
eccentricity ratio for different Re, and figure 6 depicts the same coefficient as 
a function of the Reynolds number for different orbit radius e. From the figures 
it is evident that turbulence has a large effect on the damping coefficient, and 
consequently on the tangential force. This is due to the increase in the apparent 
viscocity of the fluid as the inertia parameter grows. 

Figure 7 shows the direct inertia coefficient D /Re for various orbit 
radius, the pattern of the curves is the same as in figure 3 for moderate Reynolds 
numbers. It is clearly seen that the effect of turbulence is to increase the 
coefficient, especially at large whirling orbits, this is due to the increase in 
flow rate q fl as e grows in order to satisfy continuity of the pressure field. 

__ _A comparison of the results given in figures 5 and 7 shows that the ratio 
D /C l ess 1/10 for all eccentricities and Reynolds numbers considered, 

tHs may be an important result since it implies that the tangential force will be 
larger than the purely inertial radial force. 


Short Journal Bearing Assumption 


For the short SFD assumption, the axial pressure gradient equation comes to be: 


|| - -(12+0. 005^ (L/b) lqel) 
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H 


rU 

L 10 30 


44 3_ 2 1 

35H 3g q § J 


3q 

95 


= - Y 


(26) 

(27) 


Assuming that there is no high axial flow externally induced into the SFD, the 
flow q remains unchanged from the inertialess solution. Equations (26) and (27) 
are amenable of closed form integration, the details of the same are omitted for 
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brevity; it is found that the direct inertia coefficient D is the same as given in 
equation (2) of Table 2, rr 

Figure 8 shows the direct damping coefficient C as a function of the eccentri- 
city ratio for different values of the squeeze Reynolds number Re, and figure 9 
shows the same coefficient as a function of Re for various values of the orbit 
radius e. All calculations were made for an uncavitated SFD with a clearance ratio 
of 6=0.001, and a L/D ratio equal to 0.25. From the figures a significant influence 
of turbulence on the damping coefficient is seen. As expected, the larger the 
Reynolds number, the larger the dynamic coefficient and consequently the tangential 
force increases proportionally. Once again, this effect is due to the increase in 
the apparent viscocity of the fluid as the inertia parameter grows. 


SUMMARY 


The present paper has considered the influence of inertia and turbulence on the 
flow in the annular region between a whirling damper journal, describing circular 
centered orbits, and its bearing. After an analysis of the fluid-flow equation for 
the problem, the usual assumptions considering the length of the SFD are made to 
obtain the classical long and short journal bearing approximations. 

The region of flow was assumed to be continuous, i.e. no cavitation was 
allowed in the fluid. This allowed a clear analysis of the effect of inertia and 
turbulence on the fluid film forces and the dynamic coefficients. The laminar 
solution showed the importance of the inclusion of convective inertia terms in the 
equations of motion. The resulting reversal of the "added mass effect” makes it 
clear that numerical or analytical approaches that calculate the dynamic coefficients 
for large motion amplitudes in base to small perturbations about an equilibrium 
point may be largely in error. 

In the absence of empirical coefficients for the turbulent motion in squeezing 
flows, a friction coefficient based on the Poiseuille analysis of Hirs and Elrod and 
Ng was used. This may be modified by experimental results in the future, since 
the transition from laminar to turbulent motion may turn out to appear at larger 
Reynolds numbers than here considered. As suggested in reference (4), the values 
here presented should be considered as upper bounds for the actual dynamic coeffi- 
cients and as qualitative indicators of the influence of turbulence on the flow. 

If these bounds are even approached by the real case, turbulence will be found to 
have a large effect on the direct damping coefficient for squeeze film dampers. 

The present analysis should prove to be stepping stones for future developments 
that will consider SFD T s of finite extent and also the influence of inertia in the 
boundary conditions of the flow. 
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APPENDIX A 


Coordinate System and Equations of Motion 

Consider two circular cylinders of radii a and b (>a) and assume that the 
center of the smaller cylinder rotates with constant angular velocity m Lin a circle 
of radius e about the center of the larger one* The condition that the cylinders 
do not touch is 


0<e<l 

where '6 = «1 


e m 


(b-a) 


(A* 1 ) 
(A. 2) 


Here 6 and e are the clearance and eccentricity ratios* respectively* 


The first characteristic of the geometry of lubricant films that permits 
simplicication of the problem is that the thickness of the lubricant film, h, is 
very small compared to its length or to its radius of curvature. As consequences 
of this the following assumptions are made (8). 

1) The effects of the curvature of the film are negligible. 

2) The variation of the pressure across the film is small and may be neglected* 

3) The rate of change of any velocity component along the film is small when 
compared to the rate of change of this same velocity component across the 
film and can be neglected. 

In accordance with assumption (1), we can prescribe a fixed orthogonal 
cartesian coordinate frame {x*} . in the plane of the lubricant film. See 
figure 1, where the y* axis is in the direction of the minimum film dimension. A 
moving orthogonal coordinate frame {x_^} translating with velocity T=ba) with 

respect to {x*} and its x axis perpendicular to the line joining the centers 

of both cylinders is introduced, an^ it can be shown that the flow will be steady 
to an observer moving with the , frame. 


The following dimensionless coordinates are introduced: 


e = 


X 


5-f • 


n 




T “ tCO 


(A. 3) 


where h = b6H(0), H(0) = 1 + ecos0 (A. 4) 

is the lubricant film thickness at location 0. 

Dimensionless velocity components in the two coordinate frames are defined as 


u 


u* — . v — w* — — u 

k = yk — = -w- \fjk =s = w. u = 

b u * 6b a) * Lu ’ b(u 


(A. 5) 


u* = u + 1 


The pressure and the shear stresses are made dimensionless according to 
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>-p , a .) s 2 k ~ L T = Tx y _ = Iel (v/t) 

ojv/p 5 0n wp/fi * £n wy/d ^ ■ 


(A. 6) 


and k=l for long bearing assumption 
2 

k = (L/b) for short bearing assumption 


(A. 7) 


With these considerations, the momentum and continuity equations for the 
tridimensional flow, expressed in the moving coordinate system are 


9u 2 . 1 a , — „ -2. . YU 2 . 3 
30~ + H 3^ [UV_T1YU ] + ~ + H (UW) 


= _ k i£. + I i_ x 

K 30 h 3n 9n 


(A. 8) 


3w . 1 3 


isd 4. 


+ £ tr [wv-nywu] + (uw) - -k cf) ^ + 


3? H 3n 


H 30 




l 3? h 3n 


(A. 9) 


!f + H^ [v_nYu] + H i+ l| = o (A - io) 

where y = 3H/30 and 

Re = — b 2 (A. 11) 

V 

is the squeeze Reynolds number. 

The boundary conditions appropriate for the flow are 

at n=0 u=-l,v=^=0 ( A .12) 

n=l u=-l, v=“Y, w “0 

(Note that we have neglected the velocity component due to the motion of the 
surface n=l in the 9„direction since it is of order 6). 

The pressure must satisfy appropriate conditions at the ends of the SFD and 
must be single-valued and periodic in the circumferential direction, i.e. 

4 H de =0 (A. 13) 

Equations (A. 8) to (A. 10) are integrated across the film to obtain 


Re )39 100 + 35 I 65 


H + 4T en 


Ke I? *95 + I? hE * (b/L)2 It + 4t Sri 


3 , 3 

30 q 0 + 3? q 5 


where q_ and q are the dimensionless local flow rates in the 0 and £ directions, 
“ £ 
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(A. 17a) 



^ For the type of flow considered, the range of variation of the coefficients 
{ a i j i_i is between 1.2 and 1.0 for small Reynolds numbers and large Reynolds 
numbers, respectively; thus it may be assumed that averaged values will suffice to 
obtain meaningful results. 
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Direct inertia coefficient, 



Eccentricity, orbit radius 

Figure 3. — Direct inertia coefficient D rr /Re for circular centered orbits 
as function of eccentricity — long squeeze film damper assumption; laminar 
flow solution. 
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Direct inertia coefficient 



Eccentricity, orbit radius 

Figure 4. - Direct inertia coefficient D rr /Re for circular centered orbits 
as function of eccentricity - short squeeze film damper assumption; laminar 
flow solution* 



Direct damping coefficient. 


IP 


10 


10 

Eccentricity, orbit radius 

Figure 5. - Direct damping coefficient D tt for circular centered orbits as 
function of eccentricity - long squeeze film damper assumption; turbulent 
motion solution. 



Direct damping coefficient, 



10 10 10 10 

Reynolds number. Re 

Figure 6. - Direct damping coefficient C tt for circular centered orbit as 

function of Reynolds number - long squeeze film damper assumption; turbulent 
motion solution. 





Direct damping coefficient. 



Eccentricity, orbit radius 

Figure 8. - Direct damping coefficient c tt for circular centered orbits as 
function of eccentricity - short squeeze film damper assumption; 1/0=0.25; 
turbulent motion solution. 
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Direct damping coefficient, 



Reynolds number. Re 

Figure 9. - Direct damping coefficient c tt for circular centered orbits as 
function of Reynolds number - short squeeze film damper assumption; 
L/D=0.25; turbulent motion solution. 
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